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ABSTRACT 


A technique  is  described  for  determining  D-region  electron  density  profiles 
from  VLF  reflection  coefficients.  Some  of  the  problems  concerning  convergence 
of  the  iterative  scheme  are  discussed  as  well  as  the  constraints  introduced 
to  account  for  features  of  the  profile  for  which  information  is  not  provided 
by  reflection  coefficient  data.  Also,  an  example  is  given  where  the  technique 
is  applied  using  simulated  data. 
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I.  INTRODUCTION 

Radio  sounding  at  very  low  frequencies  (VLF  - 3 to  30  kHz)  is  a valuable 
method  for  exploring  the  lowest  ionosphere.  With  special  sounding  techniques, 
such  as  described  in  reference  1,  it  is  possible  to  obtain  continuous  sound- 
ing data  simultaneously  at  many  frequencies  over  the  VLF  band.  However,  con- 
siderable difficulty  is  encountered  in  attempts  to  obtain  quantitative 
information  concerning  ionospheric  properties  from  the  sounding  observations. 
The  theory  of  radio  wave  reflection  from  an  inhomogeneous,  anisotropic  plasma 
such  as  the  ionosphere  is  rather  well  developed,  so  that  if  the  properties  of 
the  ionosphere  are  known,  it  is  possible  to  calculate  the  expected  effects  on 
the  radio  waves.  The  real  problem  to  be  solved,  however,  is  the  inverse 
problem  of  determining  the  unknown  ionospheric  properties  from  the  observed 
effects  on  the  radio  waves.  The  use  of  data  in  this  manner  to  obtain  electron- 
density  distributions  of  the  ionosphere  has  become  known  as  profile  inversion. 

Early  work  in  ionospheric  inversion  (e.g.  references  2 and  3)  used  trial 
and  error  techniques  to  deduce  profiles.  That  is,  values  for  the  ionospheric 
properties  were  assumed  and  the  ionospheric  reflection  coefficients  calculated; 
these  were  then  compared  with  experimental  values,  and  if  the  agreement  was 
not  satisfactory,  the  calculations  were  performed  again  with  different  values 
for  the  ionospheric  properties.  The  goodness  of  fit  between  the  calculations 
and  the  experimental  data  were  determined  only  by  a subjective  estimate  of  the 
Investigator. 

There  has  long  been  a question  as  to  whether  VLF  radio  propagation  data 
might  be  readily  inverted  to  obtain  D-region  electron  density  distributions, 
without  the  use  of  trial -and-error.  There  has  also  been  a question  as  to 
whether  profiles,  if  found  to  fit  the  data,  could  be  claimed  to  be  unique. 

In  this  report  a method  of  search  is  proposed  for  investigating  the  problem. 
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In  the  case  of  steep-incidence  sounding  the  measured  data  are 
elements  of  the  ionosphere  reflection  matrix  at  several  experimental 
frequencies.  If  an  electron  density  profile  is  known,  or  assumed,  the 
elements  of  the  reflection  matrix  may  be  computed  in  a straightforward 
manner  by  use  of  a full-wave  solution,  such  as  that  given  by  Budden, 
reference  4.  The  inversion  problem  is  to  begin  with  reflection  coeffi- 
cient data  and  deduce  the  electron  density  profile. 

Many  of  the  fundamental  ideas  of  the  inversion  technique  described 
in  this  report  were  previously  presented  in  reference  5.  For  this 
reason,  several  references  will  be  made  to  that  report  rather  than  to 
reproduce  the  material. 

The  approach  taken  in  searching  for  a best-fit  electron  density 
profile  from  data  is  that  the  profile  is  to  be  linear  (on  a log  scale) 
unless  the  data  gives  information  to  justify  a given  degree  of  detail. 
To  carry  out  this  approach,  two  functions  are  defined.  The  first  is  a 
measure  of  deviation  of  computed  reflection  coefficients  from  data, 
defined  as: 
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m 

E 


where 

m - is  the  number  of  propagation  frequencies 

r^  - are  the  computed  values  of  the  reflection  coefficients 

- are  the  reflection  coefficient  data 

- are  the  uncertainties  in  the  data  values. 
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The  second  function  is  a measure  of  detail,  or  curvature,  in  the  profile  and 
is  defined  as: 


C - 


/ 


d log  Ng 


dz 


dz 


where  N.  is  the  electron  density,  z is  the  height  variable  and  the  integral 
is  taken  over  the  range  of  the  profile. 

The  required  profile  is  sought  by  an  iterative  technique  in  which  the 
total  curvature  in  the  profile  is  allowed  to  increase  in  steps  to  improve 
the  fit  to  data.  Each  step  requires  a perturbation-type  solution  in  which 
the  full-wave  solution  of  Budden  is  analytically  differentiated  with  respect 
to  perturbations  in  the  profile. 

The  end  result  sought  in  using  the  inversion  procedure  is  the  deter- 
mination of  a profile  which  contains  all  the  detail  justified  by  the  data, 
but  which  does  not  contain  spurious  detail  which  would  represent  fitting  to 
error  in  the  data.  There  should  be  an  optimum  trade-off  between  deviation 
from  data,  represented  by  the  value  of  s,  and  curvature  in  the  profile, 
represented  by  the  value  of  c.  The  trade-off  may  be  represented  by  the 
condition: 

c + Xj“nnin  (d 


where  the  value  of  X must  be  chosen  for  the  optimum  condition. 

A solution  to  (3)  may  easily  be  found  for  X»0.  A solution  which  Is 
first  order  in  r^  may  then  be  used  to  obtain  a solution  for  a small  value 


of  X.  The  first-order  solution  may  then  be  used  to  find  profiles  for 
successively  larger  values  of  X.  As  the  search  proceeds,  the  profile 
will,  in  general,  develop  more  detail  and  the  fit  to  data  will  improve 
until  a singularity  is  encountered  in  the  solution. 

This  report  is  divided  into  three  major  sections  other  than  the 
Introduction.  Section  II  is  a discussion  of  the  general  background 
needed  for  understanding  the  ideas  of  smoothing  and  inversion.  Section 
III  describes  the  procedure  for  obtaining  the  electron  density  profile 
from  reflection  coefficient  data.  Section  IV  summarizes  the  inversion 
scheme  and  gives  an  example  using  simulated  data. 

The  notation  used  throughout  this  report  is  as  described  below.  In 
general  the  subscript  convention  used  is  that  j and  j'  refer  to  the  "n" 
layers  of  the  electron  density  profile  while  i and  i' refer  to  the  "m" 
data  parameters.  Note  that  most  equations  are  written  so  that  subscripted 
variables  may  be  thought  of  as  matrices  or  vectors. 

The  notation  is: 
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Column  vector: 


II.  GENERAL  BACKGROUND 

A.  Inversion  Without  Smoothing 

A simple  Inversion  technique  Is  Illustrated  by  the  Newton-Raphson 
Iteration  procedure.  The  problem  Is  to  find  "x"  such  that: 


■ffx)  =:  r (i.e.  a constant:) 


Consider  figure  1 as  an  example  where  the  curve  for  f(x)  Is  given 
by  the  arc  PQ.  The  solution  for  "x"  Is  obtained  by  starting  with  an 
initial  value  for  "x"  (e.g.  x = x^).  This  gives,  from  figure  1: 


F-i(Xo)  - di 


rX,-  X.) 


c/Xlx=  X, 


(O 


or  let 


= (F -icx,)) 


C7) 


and 


This  leads  to  an  approximate  value  of  "x"  for  f{x)  « F,  given  by; 

= X,+  (?) 
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The  Iterative  procedure  is  then  repeated  with: 


and 


if 


(“) 


until  convergence  occurs  when  Af  becomes  less  than  some  pre-assigned 
tolerance  value.  Figure  1 illustrates  the  iterative  process. 

The  technique  will  yield  the  wanted  solution  for  "x"  provided 
that  the  slope  of  the  curve  does  not  become  zero  along  the  arc  QP  of 
figure  1. 

Figure  2 shows  an  example  of  a curve,  f(x),  where  the  inversion 
technique  will  not  converge  because  of  the  zero  slope  of  the  curve 
at  the  point  x^. 

Next,  a modified  version  of  the  Newton-Raphson  iteration  scheme  is 
considered.  In  this  case  it  Is  desired  to  follow  the  curve  f(x),  from 
an  initial  chosen  value  of  "x",  as  f(x)  approaches  the  fixed  value,  F. 
This  situation  is  Important  in  that  it  serves  as  an  analogy  for  more 
complex  cases. 
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The  problem  is  to  find  "x"  such  that  f{x)  = F (a  constant).  The 


solution  for  the  i-th  step  is: 


^ ] - 1 


where 

i(>;)  = F - ( F - ■ e" 


and  where  "X^"  takes  on  values  monotonically  increasing  from  0 to 
It  is  seen  from  figure  3 that  the  steps  converge  to  the  value 
f(x)  = F only  if  f(x)  is  a monotonic  function  of  "x".  If  this  is  not 
the  case,  the  steps  will  not  converge.  Even  in  this  case,  however, 
there  may  exist  another  function  g(x)  = function  (f(x))  which  is  a mono- 
tonic function  of  "x".  The  above  curve-following  procedure  may  then  be 
used  with  g(x)  to  converge  to  a value  G = function  (F).  This  point  is 
important  In  making  analogies  with  more  complex  inversion  problems. 

B.  Smoothing  With  No  Inversion 
1.  General  Case. 

When  data  of  the  form  y(x)  vs.  x is  considered,  it  Is  common 
practice  to  "jinooth  out"  experimental  error  by  drawing  a smooth  curve 
close  to,  but  not  necessarily  through,  the  data  points.  The  object  is 
to  draw  a curve  which  contains  all  the  detail  justified  by  the  data,  but 
which  does  not  contain  spurious  detail  which  would  represent  fitting  to 
error  in  the  data. 


Figure  4 Illustrates  a set  of  data  points  to  be  examined.  The 
method  of  smoothing,  presented  here.  Is  a technique  for  drawing  a smooth 
curve  y(x)  close  to  a set  of  data  points  = Y(X^).  In  this  technique 
the  condition  specifying  the  smooth  curve  y(x)  Is  given  by: 

c -h  >iS  — Xnin 

where  "c"  1s  the  curvature  and  "s"  Is  a measure  of  deviation  of  the 
curve,  y(x)  from  the  data  points.  The  parameter  "X"  specifies  the  trade- 
off between  fit  to  data  and  smoothness  of  the  curve  y(x). 

The  measure  of  deviation  of  the  curve,  y(x)  from  the  data  points  Is 
given  by: 

m , 

.-2-  jj  a,)  - >'a,ij-  Q.'f) 


The  Y^  = Y(X^)  are  data  values  and  the  y(X^)  are  values  of  the  curve 
y(x)  at  the  coordinates,  x = X^.  The  are  the  uncertainties  In  the 
data  values. 

In  matrix  notation,  equation  (15)  Is  written  as: 


i 


where  both  right  hand  factors  are  vectors.  Also,  1 » l,...,m;  where 
"m”  Is  the  number  of  data  points. 
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Differentiation  of  equation  (16)  with  respect  to  y(X^)  gives; 


I 


- 


\ y(x^)-  yrx-) 


.1 


(17) 


The  curvature  term  "c"  is  defined  by: 


■m 


dx 


(1?) 


Equation  (14)  defines  the  entire  curve,  y(x).  However,  in  practice, 
y(x)  can  only  be  determined  at  a finite  number  of  points,  Xj.  For 
simplicity  it  will  be  assumed  that  the  x.'s  are  equally  spaced  at  inter- 

J 

vals.  Ax,  and  that  each  X-  is  coincident  with  an  x..  From  appendix  A 

■ w 

the  second  derivative  at  x.  is  approximated  by; 

J 


__  l^h  1 ■ >>  >> ' -^7  - 1 W 

\ Ax  Ax  //  ^ 


The  curvature  function  "c"  Is  then; 


^ ’ ' 
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where  n 1s  the  number  of  points,  Xj.  This  can  be  written  in  matrix 
notation  as: 


where  the  are  elements  of  the  matrix 


1 -2 

1 

0 \ 

-2  5 

-4 

1 

\ 

1 -4 

6 

-4 

‘ ) 

\ 

\ 

\ 

\ \ ’ 

\ 

\ 

, \ 

\ \ 

\ 

\ 

\ \ \ 

\ 

\ 

\ 

\ \ \ \ 

1 

4 6-41 

0 

14  5 -2  / 

1-21/ 

The  derivation  of  this  C matrix  is  given  in  appendix  A. 

The  derivative  of  the  curvature  term,  "c",  with  respect  to  each 
y.  is  derived  in  appendix  B to  be: 


The  condition  for  the  minimization  of  equation  (14)  is: 


4. 


1 
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where  (a/3yCXj))  implies  differentiation  with  respect  to  y at  every 
point  Xj. 

From  equation  Cl7)  the  second  term  of  equation  (24)  is: 


if 

j y(  X,) 


= X, 


where 


ik)  = C^O 


y,=  x. 


Substituting  equations  (23)  and  (26)  into  equation  (24)  gives: 


where  i>l....m  refers  to  the  data  values. 


J&zV  f.A  _ 


= P (aO 


while  j»l,...n  refers  to  the  points  yj  (or  y(Xj))  which  describe  the 
the  smooth  curve  y(x). 


Next,  various  other  values  are  assigned  to  the  "X"  parameter  and 
equation  (27)  Is  solved  for  other  sets  of  y(Xj)'s.  For  any  given  set 
of  y(xj)'s  a measure  of  the  deviation  from  the  data  points  Y(X^)  can  be 
computed  from  the  equation: 


yq,)- ygj . 


y(x)  - YU-y 


r 


Therefore,  the  set  of  points  yCx.),  as  determined  for  a given  value  of 

sJ 

the  "X"  parameter,  will  be  characterized  by  a value  of  the  deviation 
term,  "s",  as  computed  by  equation  (28). 

The  following  example  illustrates  the  smoothing  technique.  Consider 
measurements  of  the  electron  density  of  the  Ionosphere  vs.  height  as 
data  values  obtained  from  a vertically  Incident  rocket.  The  problem  Is 
to  "smooth  out"  the  error  in  the  data  and  interpolate  between  data  points 
by  drawing  a smooth  curve  close  to  but  not  necessarily  through  the  data 
points.  The  error  in  each  measurement  will  be  assumed  to  be  gaussian  with 
standard  error,  o.. 

There  are  20  heights  at  which  data  points  are  taken  and  these  height 
values  are  2 km  apart.  The  height  values  for  the  computed  smooth  curve 
will  be  1 km  apart.  To  apply  the  procedure  as  presented  by  equat1ons(14) 
through  (28),  let  the  height  parameter  be  taken  as  the  Independent 
variable,  x,  and  the  log  of  the  electron  densities  (I.e.  log^Q  N)  be 
taken  as  t;ie  dependent  variable,  y.  The  optimum  value  of  the  deviation 
parameter  "s"  Is  obtained  when  the  fit  between  the  data  and  the  curve 


14 


y(x)  = 1o9io  average,  within  one  "a"  value  at  each  data 

point.  That  is,  for  20  data  points  "s(optimum)"  will  be  equal  to  20. 

Figures  5 through  9 illustrate  the  smoothing  procedure  as  applied  to 
the  rocket  experiment.  Figure  5 shows  what  will  be  called  the  true 
electron  density  profile  along  with  the  "data"  values  which  are  obtained 
by  adding  random  gaussian  values  to  the  true  density  values  at  each  of 
the  20  heights. 

Figure  6 presents  the  results  obtained  after  applying  the  smoothing 
technique  to  the  data  (with  gaussian  error  included).  Shown  are  the 
resulting  curve  for  "optimum  smoothing"  (i.e.  s=r20)  and  the  curve  for 
"over  smoothing"  (i.e.  s = 31.1).  Over-smoothing  implies  that  minimiza- 
tion of  the  curvature  term,  "c"  of  equation  (14)  is  being  accomplished 
at  the  expense  of  not  fitting  well  to  the  data  points.  That  is,  the  choice 
of  A is  too  small  and  needs  to  be  increased. 

Figure  7 compares  the  true  profile  with  the  optimumally  smoothed  profile. 

Figure  8 gives  curves  for  "optimum  smoothing"  and  for  the  case  of 
"under-smoothing"  (i.e.  s = 9.1).  In  the  instance  of  "under-smoothing" 
the  minimization  of  the  term  for  deviation  to  the  data  points  (i.e.  "s") 
is  dominant  over  minimizing  the  curvature  term.  The  smooth  curve  y(x)  is 
being  attracted  more  strongly  to  the  data  points.  This  result  implies 
that  the  choice  of  "X"  is  too  large  and  should  be  decreased. 

If  the  value  of  "X"  is  increased  still  further  than  that  of  figure  8, 
it  is  shown  in  figure  9 that  the  resulting  smooth  curve,  y(x),  is  a fit 
to  the  "error"  rather  than  to  the  true  profile.  This  result  is  an 
important  aspect  of  the  smoothing  procedure  when  there  is  error  in  the 
data  values. 


% 
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2.  Limiting  Case  of  Interpolation  Between  Data  Points. 

Given  a set  of  data  points  Y(X^),  where  1«l,...m.  Next  consider 
the  set  of  evenly  spaced  points  Xj,  j“l,...n,  which  Include  the  set  X^. 
The  problem  Is  to  find  the  values  of  terms  y(Xj)t  corresponding  to  the 
various  Xj's,  by  interpolation  between  the  Y(X^)  values. 

Consider  equation  (27),  that  Is: 


fAxr"  i 


l\T 


07"  V 


. o (XT) 


In  the  limit  as  "X" -*«»  equation  C27)  becomes: 


J , 

JA / j\  / JA 


(a<f) 


where: 


yj  » are  interpolated  values  at  evenly  spaced 

Increments  In  x (I.e.  and  which  Include  the 
set  of  Y(X^)'s. 

Y^  * Y(X^),  1 ■ l,...m  are  data  values  at  each  X^. 

Note  that  each  Y^  must  be  Identical  to  a yj.  Also  the  modified  C matrix 
Is  Identified  by  the  superscript  Il-b  to  differentiate  It  from  other 
forms  used  in  later  sections  of  this  report. 

Il-b 

The  J'-th  row  of  the  modified  matrix,  Cjij  , Is  the  j'-th  row  of 

the  matrix  C,,i  (see  equation  (22))  If  Xj,  does  not  correspond  to  an  X^. 

j J JJ.5  ' 

If  Xj)  does  correspond  to  an  X^,  then  the  j'-th  row  Is  all 
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Also,  a11 


zeros  except  for  a 1 for  the  main  diagonal  element,  Cjiji. 
elements  of  the  vector  Vj,  of  equation  (29)  are  zero  except  that  If 
Xj,  corresponds  to  an  X^,  then  Vj,  = Y^. 

As  an  example,  consider  the  set  of  points  Yj,  Y2  and  Y^  at  Xj,  X2 
and  Xj.  Now  assume  there  are  a set  of  Xj's  (separated  by  even  Incre- 
ments) positioned  between  the  X^'s.  Let  there  be  4 Xj  Increments 
located  between  successive  X^  values  with  two  Xj  Increments  preceding 
Xj  and  two  increments  of  Xj  following  Xj.  That  Is,  x^  = Xj,  x-j  = X2, 
and  Xjj  = Xj.  When  these  assumptions  are  substituted  Into  equation  (29), 
the  following  matrix  expression  results. 


o 
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The  Interpolated  values,  y^,  are  then  found  from  the  matrix 
equation: 


It  Is  shown  In  appendix  C that  the  Yj's,  obtained  by  using  this 
smoothing  technique,  are  linear  combinations  of  the  Y^'s.  That  Is: 


It  is  also  shown  In  appendix  C that  derivatives  at  the  interpolation 


points,  X.,  are  given  by: 

V 


m 


s.  --  L y; 


(i!-) 


h-1 


where 


00 

So  that  the  curve,  r(w)  is  a function  of  the  curve,  a(z). 

Now,  suppose  that  at  a finite  number  of  values  of  w (i.e. 

1=1,... m)  there  are  data  values,  R(Wj)  = value^*^!^  * ^i*  figure  10. 

The  error,  is  unknown;  however,  the  expected  error  is  a^.  From  these 
data  values,  R(w^),  it  is  desired  to  deduce  the  form  of  the  curve,  a:(z), 
insofar  as  possible.  The  curve  a(z)  can  not  be  determined  exactly 
because  (a)  there  are  only  a finite  number  of  data  values,  R(cj.)  and 
(b)  there  is  unknown  error  in  the  data.  This  problem  is  one  of  both 
smoothing  and  inversion. 
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C.  Combination  Problem  of  Smoothing  and  Inversion 
Define  a function  r{co)  as: 

b 

f(u))  — j |^K(  u;  o<r/ij 


Because  of  the  error  in  the  data,  R(<*>j),  it  may  not  be  desirable 
to  represent  a{z]  as  a curve  for  which  the  r(Wj)  fit  the  data  values 
exactly.  The  optimum  curve,  a(z),  should  be  one  which  contains  all  the 
detail  Justified  by  the  data,  but  which  does  not  contain  spurious 
detail  which  would  represent  fitting  to  error  in  the  data. 

The  trade-off  between  fit  to  data  and  smoothness  of  the  curve  may 
be  represented  by  equation  (14)  as: 

C -h  5 — ^ /7?/t7 /rT^n/yj  Q.  ¥ 

where  "c"  is  the  measure  of  curvature  in  a(z)  and  "s"  is  a measure  of  the 
deviation  of  values  r(co^),  on  the  curve  r(w),  from  the  data  values. 

The  function  "s"  is  defined  as: 


where  the 

o^'s  are  the  uncertainties  in  the  data. 
The  function  "c"  is  defined  as: 
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The  parameter  "X"  of  equation  04}  is  chosen  so  as  to  optimize  the 
trade-off  between  the  fit  to  the  data  RC^^)  and  the  smoothness  of  the 
curve  “(z). 

For  the  sake  of  numerical  tractability,  the  curve  a{z)  is  repre- 
sented by  a series  of  short  segments,  see  figure  11.  The  values  of  «(z) 
where  the  segments  join  (i.e.  at  Zj,  j=l,...n)  are  taken  to  be  the 
unknown  parameters  of  a(z)  and  are  denoted  by  «(Zj)  = aj,  j=l,...n. 

Because  a(z)  can  only  be  determined  at  a finite  number  of  points, 

Zj,  and  since  for  simplicity  it  will  be  assumed  that  the  Zj's  are  equally 
spaced  at  intervals  Az,  the  derivative  in  equation  (39)  at  Zj  is 
represented  by: 


^ 2 


As  given  previously  in  equation  (20),  the  curvature  function,  "c"  can 
then  be  approximated  by: 


where  n is  the  number  of  segment  points,  Zj. 

Equation  (41)  can  be  written  in  matrix  notation  as: 
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where  the  matrix  Cj,j  is  given  1n  equation  C22). 

The  minimization  of  equation  (14)  Is  with  respect  to  the  positions 
of  the  points  “(zj)  on  the  curve  a(z).  Therefore,  the  minimization 
condition  Is  met  If: 


-t-  A ii 


where  j=l,...n 

The  derivative  (3c/3a(z.))  Is  given  by: 

J 


which  Is  similar  to  equation  (23). 

The  derivative  (ds/3a(zj)  Is  given  by: 

J 


where  r^  » r(u;^),  (1=1. ..m)  Is  computed  from  equation  (37)  using  the 
values  ttj  » a(Zj).  For  purposes  of  Integration,  linear  Interpolation  Is 
made  between  the  points  "(Zj)> 
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The  derivative,  (ar^/aaj),  Is  given  by: 


=/  { 


4- 


j« 


(4^ 


where  (3a(z)/9aCzj))  Is  a triangular  weighting  function  made  necessary 
by  the  use  of  linear  Interpolation  in  the  Integration  and  is  given  by: 

4 


,)o(Cz)  — I 


-Ji-i  \ '•*■  Ii-1  < ? i [(^7) 


O^her  aJise^ 


See  figure  12,  where  the  ratio  (3a/3a.)  is: 


— — ~ \ 


Substituting  equations  (44)  and  (45)  into  equation  (43)  gives: 


where  r^  = r(<«>^)  and  Is  given  by  equation  (37).  Also,  the  rectangular 
matrix  (ar^/^a^^  Is  given  by  equations  (46)  and  (47). 
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Note  that  the  unknowns  in  equation  (48)  are  the  aj's  but  that  the 
r/s  are  functions  of  these  "unknown"  « 's.  Hence,  at  this  point 
equation  (48)  cannot  be  used  to  solve  for  the  a(z)  curve. 

With  certain  additional  constraints  in  the  curve,  a(z),  the  equation 
may,  however,  be  solved  for  "X"  = 0,  for  which  values  of  r{u>.)  are  not 
required.  Without  additional  constraints  the  solution  for  "X"  = 0 is 
indeterminate  since  any  straight  line  form  of  a(z)  is  a solution.  That 
is,  for  any  straight  line,  the  curvature  function  "c"  = 0.  The  required 
additional  constraints  could  be,  for  example,  that  the  values  a(a)  = 
and  a(b)  = be  pre-selected  fixed  values.  These  constraints  may  be 
added  by  modifying  certain  elements  of  the  C matrix  of  equation  (22). 

That  is  by  setting: 


c; 


h h-i 


o 


and  by  setting  the  first  element  of  the  right-hand  vector  of  equation  (48) 
to  Q(a)  and  the  last  element  toa(b).  If  these  constraints  are  retained 
for  other  values  of  "X"  as  well,  equation  (48)  becomes: 


! 

i 

i 
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II-C 

where  the  modified  C matrix,  C contains  the  changes  just  described  and 
contains  the  factor  (1/Ax  ).  Also,  for  this  case  the  vector  "V"  is 


Note  that  other  constraints  are,  of  course,  possible. 

The  solution  to  equation  (50)  for  "X"  = 0 is  a set  of  aj's  repre- 
senting a straight  line  between  the  fixed  points  a(a)  and  «(b).  Although 
exact  solutions  to  equation  (50)  may  not  be  obtained  for  other  values 
of  "X",  a first-order  solution  may  be  found  for  a small  positive  value 
of  "X"  if  r^  is  approximated  by  a first  order  Taylor  expression: 


where  a refers  to  the  solution  for  X = 0 and  r^°  and  (ar^/aaj  are 
values  obtained  from  equations  (37)  and  (46)  respectively.  Equation  (50) 
is  then  aporoximated  by: 
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If  A Wj  is  defined  as  C“j  - “j°)»  equation  (53)  may  be  written 


In  this  matrix  equation  only  the  Aa.'sare  not  known  and  hence  the  equation 

V 

may  be  used  to  solve  for  Aa  Then: 


— o(  C ^j)  +•  j I r h 

is  an  approximate  solution  for  the  small  value  of  "X"  used  in  the  equation. 

If  the  ttj®  values  in  equation  (54)  are  now  taken  to  be  the  values 
just  found  for  small  "X",  this  equation  may  be  used  with  a somewhat  larger 
value  of  "X"  to  obtain  a new  curve,  a(z),  represented  by  values,  a..  Thus, 

s) 

as  long  as  the  increments  in  "X"  are  sufficiently  small  and  so  long  as  tl.e 
coefficient  matrix  does  not  become  singular,  equation  (54)  may  be  used  in 
"bootstrap"  fashion  to  obtain  values  of  Oj  for  any  value  of  "X".  That 
is,  for  a given  value  of  "X",  equation  (64)  is  solved  for  (Aa^).  Then 
(new)  is  computed  from  equation  (55)  as  Oj  (new)  Qj  (old) 

At  this  point  new  values  of  r^  and  (dr^/d«j)  are  computed.  The  next  step 
is  to  substitute  these  "new"  values  Into  equation  (54),  along  with  a new 
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value  of  "X"  and  then  to  solve  the  equation  for  yet  another  set  of 

(A.j). 

Equation  (54)  may  also  be  written  In  a simplified  form  as: 

(Ki  • - m (^0 

where  [K]  Is  fte  coefficient  matrix, 

(Aotj)  Is  the  unknown  vector 
and  (T)  Is  a known  vector. 

The  solution  for  (Aa.)  is  found  from: 

J 


It  may  be  noted  that  equation  (57)  Is  similar  In  form  to  equation  (12). 

In  section  A It  was  shown  that  the  modified  Newton-Raphson  procedure 
diverged  If  the  derivative  of  f(x)  became  zero  anywhere  along  the  f(x) 
curve  being  followed.  That  Is  where: 


An  analogy  exists  between  equation  (57)  and  equation  (58)  In  that 
the  procedure  for  solving  equation  (57)  will  also  diverge  If  the  coeffi- 
cient matrix,  K,  becomes  singular. 

Also,  In  analogy  with  section  A,  If  divergence  occurs,  there  may  be 
another  function,  g(uj)  = function  (r(cj)),  which  may  be  followed  to  the 
desired  Inverse  solution  without  a singularity  being  encountered.  The 
choice  of  function,  g(w),  however,  is  likely  to  be  an  art  more  than  a 
science,  unless  sufficient  information  is  available  on  the  behavior  of 
these  functions. 

The  matrix  equation  to  be  used  with  g{w)  = function  (r(<j))  Is 
simply  equation  (54)  with  the  r^'s  replaced  by  g^'s  and  the  appropriate 
choice  made  for  the  matrix  of  uncertainties,  o^.  These  "uncertainties" 
are  probably  best  determined  by  requiring  that  as  "X"  Is  Increased  to 
make  the  values  of  g^  approach  the  values  = function  (R^),  the  sequence 
of  curves,  a(z),  become  the  same  as  If  the  function  r(uj)  were  used.  The 
formulation  needed  for  making  such  a choice  of  "uncertainties"  in  G(cOj) 

Is  given  In  appendix  D,  where  it  Is  seen  that  the  uncertainty  matrix  for 
the  G^'s  Is  no  longer  diagonal.  For  the  g^'s  equation  (54)  becomes: 


I 
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III.  DETERMINATION  OF  ELECTRON  DENSITY  PROFILES  FROM  VLF  REFLECTION 

COEFFICIENTS 

A.  Reflection  Coefficient  Data 

Very  low  frequency  (VLF)  reflection  coefficient  data  may  be  obtained 
from  steep- incidence  sounders,  such  as  the  NOSC  facility  located  at  Sentinel, 
Arizona,  reference  6.  This  system  uses  a horizontal  dipole  as  a transmitting 
antenna  and  crossed  loops  for  receiving.  The  system  has  the  ability  to  trans- 
mit several  frequencies  nearly  simultaneously. 

Reflection  of  a VLF  wave  at  near-vertical  incidence  is  described  by  the 
reflection  matrix: 


where  the  first  subscript  refers  to  the  polarization  of  the  upgoing  wave  with 
respect  to  the  plane  of  incidence,  and  the  second  subscript  refers  to  the 
polarization  of  the  downcoming  wave. 

The  measured  parameters  of  the  sounding  system  consist  of  the  received 
electric  fields,  both  parallel  and  perpendicular  to  the  plane  of  incidence, 
of  the  VLF  radio  wave.  These  fields  are  identified  as  E,,  and  E^  respectively. 

The  transmitted  wave,  Ej,  is  used  to  form  the  ratios: 
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Equation  (65)  may  also  be  re-written  in  the  form: 


Rc  cUi) 


C^O 


where  i = 1,. . .m 

In  this  case  the  terms  are  identified  as: 

= The  transmitted  frequency. 

m = The  total  number  of  transmitted  frequencies. 

R(caj)  = The  value  of  the  ratio,  at  the  i-th  transmitted  frequency,  as 
determined  from  data. 

6 - The  incident  angle  of  the  radio  wave  on  to  the  ionosphere. 

j = \/^ 

It  is  to  be  noted  that  the  above  reflection  coefficients  are  functions  of 
the  earth's  magnetic  field.  In  particular,  as  the  location  of  the  sounder  site 
approaches  the  geomagnetic  equator,  the  cross  term  (i.e.  tends  to  zero 
and  the  ratio  term,  R(uj^)  becomes  equal  to  one.  When  this  situation  occurs 
the  profile  determination  technique  then  yields  no  information.  The  most 
useful  data  for  determining  profiles  is  obtained,  therefore,  at  locations  with 
large  magnetic  dip  angles  (i.e.  near  the  geomagnetic  poles). 

In  equation  (66)  all  of  the  reflection  coefficient  terms  imply  complex 
numbers.  Also,  the  ratio  term,  R(oj^)  is  a complex  number.  A useful  repre- 
sentation of  the  reflection  coefficient  data  is  to  plot  it  in  the  complex  R- 
plane  as  a function  of  the  data  frequencies,  An  example  of  this 
relationship  is  Illustrated  in  figure  13  for  the  real  and  imaginary  components 
of  R(ui|).  Note  that  the  frequency,  w,  serves  as  the  parametric  variable. 

Included  in  figure  13  is  a smooth  curve  which  would  represent  the  data 
if  it  were  obtained  at  all  frequencies  without  experimental  error. 


B.  The  Ionospheric  Model 

It  1s  proposed  to  investigate  the  possibility  that  the  form  of  the  elec- 
tron density  profile  of  the  ionosphere  can  be  deduced  from  steep  Incidence 
sounder  measurements  of  reflection  coefficients  (equation  (62)),  since  the 
measured  values  depend  on  the  height  variation  of  electron  density.  These 
values  also  depend  on  the  charge,  and  mass  of  the  particles  In  the  Ionosphere, 
as  well  as  their  collision  frequencies  with  neutral  particles  and  of  the 
earth's  magnetic  field. 

Several  assumptions  about  the  electron  density  vs.  height  relationship 
In  the  Ionosphere  must  be  made  before  a mathematical  technique  can  be  set  up 
to  determine  a profile  from  reflection  coefficient  data.  In  particular,  all 
of  the  above  parameters  will  be  taken  as  known  except  the  electron  density 
profile  which  will  be  considered  as  unknown  and  Is  to  be  found  from  the 
da  ta. 

The  magnetic-field  parameters  will  be  taken  to  be  known  and.  In  the  case 
of  the  model  presented  here,  the  magnetic  azimuth  will  be  restricted  to  be 
either  90*  or  270°.  As  is  the  usual  case  of  these  VLF  frequencies.  Ions  are 
considered  to  have  negligible  effect  on  propagation  and  hence  are  omitted  from 
the  model.  The  collision  frequency  profile  Is  taken  to  be  known  and  Is 
constrained  to  be  of  exponential  form  (e.g.  *>  = Exp  (-az)). 

Additional  assumptions  (i.e.  constraints)  on  the  form  of  the  desired 
electron  density  profile  are  listed  below  and  are  Illustrated  In  figure  14 
where  the  electron  density  (Ng)  Is  shown  as  a function  of  ionospheric  height 
In  terms  of  log^Q  N^. 

1.  "Stop"  constraints  are  Introduced  at  both  the  top  and  bottom  of  the 
electron  density  profile.  These  correspond  to  limiting  values  of 
electron  density  which  the  profile  may  assume  at  these  heights.  At 
the  top  of  the  profile  the  "stop"  defines  the  minimum  value  of 
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electron  density  which  the  profile  may  assume.  At  the  bottom,  the 
"stop"  defines  the  maximum  value  for  that  height. 

2.  It  is  assumed  that  there  are  no  electrons  below  the  chosen  lowest 
height  of  the  profile.  Above  the  chosen  top  height  of  the  profile 
the  electron  density  is  assumed  to  be  of  constant  value  equal  to  the 
electron  density  at  the  chosen  top  height.  This  latter  property  is 
usually  referred  to  as  a semi-infinite  medium. 

3.  The  profile  is  constrained  to  be  reasonably  smooth  but  otherwise 
consistent  with  the  reflection  coefficient  data.  Note  that  this 
implies  that  the  profile  approaches  a slope  of  (dN/dz)  = 0 near  the 
top  height  which  is  just  below  the  semi-infinite  medium. 

C.  Computation  of  Ionospheric  Reflection  Coefficients 

Ionospheric  reflection  coefficients,  R,  may  be  obtained  for  a given 
electron  density  profile  by  the  "full-wave"  method  given  in  reference  4. 
These  reflection  coefficients  are  the  result  of  the  integration: 


R- 


L 


where  the  integrand  is  given  in  reference  4 as: 


Jz  = -ik/Z  I + 


k is  the  wave  number. 

The  matrices,  "S",  are  functions  of  the  electron  density  profile  and  of  the 
wave  frequency,  co,  as  well  as  other  parameters  described  In  section  B. 

Equation  (68)  Is  given  In  scalar  form  In  reference  5,  pages  11  and  12. 
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As  was  discussed  in  section  A,  a useful  function  of  the  reflection 


coefficients  is: 


J.^11 


where  i = 1, . . . ,ni 

with  m being  the  total  number  of  data  frequencies. 

In  this  instance  the  reflection  coefficients,  j^Rj^  andj^R||  , are  obtained 
from  the  integration  in  equation  (67)  and  thus,  "r(cj.)"  is  the  computed  value 
of  the  ratio  at  the  i-th  data  frequency. 

The  function  "r(w.)"  is  complex  and  as  such  can  be  plotted  in  the  complex 
r-plane  as  a function  of  its  real  and  imaginary  parts.  Examples  of  r(cu) 
curves  are  illustrated  in  figure  15  as  a continuous  function  of  the  parametric 
variable,  cj.  Calculations  performed  using  different  electron  density  profiles 
will  produce  various  forms  of  the  r(co)  curve.  Some  of  r(co)  curves  will  possess 
certain  distinct  characteristics,  such  as  the  loop  shown  in  figure  15  for  the 
curve  |r(o;))j  ; others  may  be  relatively  smooth,  like  the  [r(w)]2  curve.  Note 
that  because  of  the  constraints  listed  in  part  B,  not  all  forms  of  the  r(cj) 
curve  can  be  generated  simply  by  varying  the  electron  density  profile. 

0.  Comparison  of  Data  with  Computed  Reflection  Coefficients 
Reflection  coefficient  data  may  also  be  plotted  in  the  complex  R-plane. 

An  example  was  presented  in  figure  13  where  the  data  points  were  plotted  para- 
metric with  frequency.  If  ionospheric  reflection  coefficients  are  computed 
from  an  assumed  electron  density  profile,  a plot  illustrating  the  comparison 
between  data  and  theoretical  calculations  may  be  constructed.  Figure  16 
shows  an  example  of  such  a comparison.  In  this  instance  it  is  observed  from 
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the  plot  that  the  computed  reflection  coefficients  differ  drastically  from 
the  data  values.  The  problem  to  be  solved,  then,  is  to  find  that  electron 
density  profile  which,  when  used  in  the  full-wave  computation  of  reflection 
coefficients,  will  produce  a match  to  the  data,  R(uj^).  Such  a match  is 
illustrated  in  figure  17.  Because  of  error  in  the  data  and  because  of  the 
constraints  put  on  the  ionospheric  model,  the  data  can  not  be  matched  exactly. 
This  characteristic  is  shown  in  figure  17. 

The  procedure  for  finding  that  particular  profile  of  electron  density 
which  will  produce  a match  to  the  data  is  discussed  in  the  following  sections 
of  this  report. 

E.  Application  of  the  Combined  Smoothing  and  Inversion  Techniques 
to  the  Ionosphere  Problem 

The  question  as  to  whether  the  the  profile  of  the  electron  density  of 
the  ionosphere  can  be  determined  from  reflection  coefficient  data  is  examined 
with  reference  to  the  procedure  presented  in  section  II,  part  C of  this  report. 

The  correspondence  between  certain  variables  may  be  identified.  In 
particular,  the  following  relations  apply: 


Section  III 

Section  II, C 

a 

Height,  z 

The  Integration  variable,  z 

b 

Wave  frequency , oj 

Parametric  variable,  w 

c 

log  Ng  vs.  height,  z 

a vs.  z 

d 

log  Ng  at  height,  Zj 

“j  = “<^4) 

e 

Re  ( r(cj)]  & Im  [ r(w)] 

r(cj) 

f 

Re  f r(w^)J  & Im  ( r(w^)J 

r^  = r(cu^) 

9 

Re  [R(ui^)]  & Im  [R(a;^)] 

R^  = R(Wj) 

Note  that  the  ionosphere  r(w)  is  not  obtained  from  an  integration  of 
exactly  the  form  given  in  equation  (37).  However,  the  integration  (67) 
together  with  equation  (69)  is  comparable  with  the  integration  in  equation  (37) 


in  that  the  r{oi)  values  are  obtained  as  a result  of  an  integration  between 
fixed  limits  in  the  integration  variable,  z,  and  furthermore  the  integrand  is 
a function  of  z,a(z),  and  w. 

The  derivative  (dr/da^)  required  in  equation  (54)  of  section  II-C  will 
correspondingly  be  somewhat  different  in  form.  It  is  first  noted  that: 


the  derivatives  (dj^Rj^/da^)  and  C3j^R||/ao(^)  are  elements  of  the  matrix 


/ ^ ll^il 

^ iRii 

' 

\ ^ iiRi 

d iRl 

\() 

which  is  found  as  the  result  of  the  integration 
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The  integrand  Is  found  by  differentiating  equation  (68)  with  respect  to 


3 


? B ' B — ^ 

d oi  ^ » 


T 


where  the  weighting  factor  (da/da.)  is  described  by  equation  (47).  The 

V 

derivatives  (3'V/3“)  are  described  in  reference  5,  pages  31-35. 

The  constraints  on  the  o!(z)  curve  when  it  is  used  to  represent  the 
electron  density  profile  are  somewhat  different  than  the  ones  placed  on  the 
a(z)  curve  in  section  II-C.  Instead  of  choosing  the  ends  of  the  a(z)  curve 
(the  profile)  to  be  fixed  values,  inequalities  are  used.  At  the  top  of  the 
profile  the  electron  density  is  constrained  to  be  at  least  as  great  as  a 
chosen  fixed  value.  At  the  bottom  height  of  the  profile,  the  electron 
density  is  constrained  to  be  no  greater  than  a chosen  fixed  value.  Thus: 

o(  Yqp  stcp 

C ^ ^ ffoT.  STOP 

In  addition,  in  order  to  ensure  a continuous  slope,  (dN^/dz),  leading 
upward  to  the  semi-infinite  medium,  the  slope  at  the  top  of  the  profile 
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(z  = a = Zj)  Is  constrained  to  approach  that  of  the  semi-infinite  medium. 
In  terms  of  a segmented  «(z)  curve,  that  is,  a segmented  electron  density 
profile,  this  is  implemented  by  adding  the  term: 


) ' iA'^) 


-3 


(75) 


to  the  summation  representing  the  curvature  function.  That  is,  "c"  is  now 
defined  as: 


= 3. 


Compare  this  result  to  equation  C20). 

With  these  modifications  to  the  set  of  constraints,  the  equation  to  be 


solved  is  a modified  form  of  equation  (54)..  In  particular  the  matrix: 

.(1  + Pp>  -1  ° ■ O 


-1 

o 


-1  o 
o 

c 0 


-f  1 


C77) 


O PC'**' 

is  added  to  the  unmodified  C matrix  of  equation  (22).  That  is,  the  original 
C matrix  terms  are  modified  to  be: 


C”!!- 

1 t ( i.  f - 

II 

1 

■4- 

1 

!-»■ 

1 1 

csi  - 

-a  t(-  1)  ; 

-3 

5-  + ( 1)  r 

io 

1 - 

ItPg  J 

00 
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Here: 


(a) 


Pt 


(b) 


Pb 


’»  0 if,  in  the  solution  for  Aaj's  the  top  end  of  the  profile 
is  Inclined  to  move  to  the  right  of  the  stop. 

= a very  large  positive  value  if  the  top  end  of  the  profile 
is  inclined  to  move  to  the  left  of  the  stop. 

= 0 if,  in  the  solution  for  the  bottom  end  of  the 

profile  is  inclined  to  move  to  the  left  of  the  stop. 

= a very  large  value  if  the  bottom  end  of  the  profile  is 
inclined  to  move  to  the  right  of  the  stop. 


V 

Also,  the  vector,  V,  of  equation  (51),  used  in  equation  (54),  is 


modified  to  be 


' r S Top 

o 


-f>  - j3oT.  sTpp 

9 A 


Substitution  of  equations  (78)  and  (79)  into  equation  (54)  gives: 


N 

(so) 


where  the  r^'s  are  the  computed  ionospheric  reflection  coefficient  ratios  as 
given  by  equation  (69),  the  R^'s  are  the  data  parameters  given  by 
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equation  (66)  and  the  oj's  are  the  logs  of  the  electron  densities  at  the 
profile  segment  points,  Zj. 

The  derivatives,  (3r^/a«j),  are  given  by  equations  (70-73). 

F . Characteristics  of  the  Solution  of  Equations  In  Terms  of  the 
"r(u;)"  Parameter 

The  solution  to  equations  of  the  form  of  equation  (80)  was  discussed 
previously  in  section  II,C  with  regard  to  equation  (54).  The  solution  for 
successive  iterative  steps  of  X+A\  is  given  in  terms  of  (Aa.).  Each 
Iterative  step  leads  to  a new  set  of  a.'s  given  by: 

Oi(Z^)„rZ  00 

The  procedure  is  to  compute  new  sets  of  the  ionospheric  parameters,  r(cj),  from 

each  set  of  a(z.) 

^ new 

In  section  D the  comparison  of  computed  values,  r(cj^),  and  data  values, 

R(to^),  were  shown  plotted  In  the  complex  r-plane.  The  tJij's  are  the  propagation 

frequencies  and  serve  as  the  parametric  variable  In  the  plots. 

It  was  suggested  that  the  proper  choice  of  electron  density  profile  (I.e. 

set  of  a^'s)  should  result  In  a set  of  r(uj(.)'s  which  would  match  the  set  of 

data  R(u;^)'s.  Figure  18  illustrates  the  iterative  procedure  corresponding 

to  successive  solutions,  [a(z)J  , to  equation  (80). 

^ new 

The  set  of  r(tj|)'s,  computed  from  each  electron  density  profile  (I.e.  each 

successive  set  of  [a(z)J  , are  presented  In  figure  18  as  the  curves 

new 

(r(cj)Ij,  [r(tj)l2»  etc.  The  data  values  are  shown  In  the  figure  as 
R(c^),  R{<^),  etc. 
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In  carrying  through  the  Iterative  procedure  to  obtain  a match  between 
the  computations  and  the  data.  It  Is  noted  the  minimization  requirement  of 
equation  (14)  implies  that  the  deviation  between  the  r and  the  R(ci:!|)'s 

be  minimized.  Following  equation  (15),  this  may  be  written  as: 

S r £ ^ 

L--1  ^ 

which  demands  that  each  frequency  point  on  the  r(tj)  curve  pursue  approximately 
the  shortest  course  possible,  in  the  complex  plane,  to  get  from  the  tr(co)]2 
curve  to  the  R(u>.)'s. 

Observation  of  the  curves,  in  figure  18,  shows  that  the  points  of 
(r(co))j  are  numbered  in  the  opposite  sense  to  those  of  the  data  R(<o^)'s.  The 
trend  of  the  r{tj)  curves  will  need  to  include  a complete  reversal  If  a match 
to  data  Is  to  be  accomplished.  The  arrows  which  emanate  from  each  r(tj^)  point 
Indicate  the  direction  in  which  each  r(cOj)  point  is  constrained  to  move  as  the 
iterative  steps  progress  toward  the  final  match  using  the  above  definition 
of  "s". 

It  Is  noted  that  in  the  sequence,  shown  in  figure  18,  that  a loop  Is 
formed  and  then  unwound.  In  forming  a loop,  however,  the  curve  Is  forced  to 
contain  a "cusp"  at  one  stage,  as  is  shown  on  the  curve  (r(w)]^.  Experience 
with  numerical  computations  has  shown  that  In  cases  of  this  type  the  coeffi- 
cient matrix  of  equation  (80)  becomes  singular  as  the  form  of  the  curve 
containing  the  "cusp"  Is  approached,  and  hence  the  solution  for  the  and 

thus  for  the  electron  density  profile,  diverges. 


G.  Transformation  Function,  "g(r,cj)" 

The  sequence  of  r(o))  curves  in  the  above  example  need  not  Include  a 
curve  containing  a "cusp"  if  one  end  of  the  r(a;)  curve  is  rotated  about  the 
other.  This  does,  however,  require  a different  definition  of  the  "s"  of 
equation  (82). 

The  definition  of  "s"  required  for  this  sequence  must  emphasize  the 
angle  assumed  by  the  curve  in  the  complex  plane  at  each  frequency  point, 
r{co. ),  except  at  the  end  about  which  the  curve  is  to  be  rotated.  At  this 
end  point  the  absolute  position  of  the  point  is  to  be  stressed,  as  before. 

A complete  definition  of  the  curve  must  also  include  a requirement  that  the 
spacing  of  frequency  points,  r(aj. ),  must  match  for  the  computed  and  data 
curves.  It  is  then  noted  that  since  the  low  frequency  end  of  the  r(tj)  curve, 
r(uj^),  is  more  stable  with  respect  to  small  changes  in  the  electron  density 
profile,  it  is  likely  to  be  the  better  choice  for  the  end  about  which  rotation 
is  to  be  carried  out. 

A function,  g{r),  containing  the  above  properties  is: 


and 


(p- 


where  4>^is  the  principle  part  of  the  phase  term  and  lies  in  the  Interval 
-w  <4>  < n . 

Figure  19  illustrates  the  transformation  function,  g,  in  the  r-plane. 
The  phase  angle,  <l>^,  may  be  written  as: 


The  value  of  n^ , in  the  phase  term,  is  chosen  for  each  frequency  so  as 
to  ensure  that  the  phase,  </>^,  is  continuous  along  the  r(to)  curve.  This 
requirement  leaves  the  value  of  n^  arbitrary  for  one  frequency,  e.g.  the 
lowest  frequency,  The  n^  value  for  this  frequency,  that  is,  n^,  is 
chosen  so  as  to  cause  the  curve  to  rotate  in  the  desired  sense  (i.e.  clock- 
wise or  counterclockwise). 

The  transformation  function  for  the  data  values,  is  introduced  as: 


4-  f R^C/7>7c,^,) 


wnere 
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The  phase  angle,  may  be  written  as: 


(4.)^  = 4RCTAW 


(87) 


Note  that  at  a "cusp",  as  in  the  curve  [ r(cj)]  ^ of  figure  18,  (dr/dw)  = 0 
and  hence: 


^ — I D ( oic4^)  ~ ^ -f  -j-i  In  ciet  b r iTtmah') 

The  new  definition  of  "s",  which  emphasizes  a small  deviation  between 
computed  values,  g.,  and  data  values,  G^,  is  given  by  equation  (61)  as: 


This  formulation  will  not  allow  a "cusp"  to  form  in  the  iteration  scheme. 

The  solution  for  the  and  hence  the  set  of  (“j)pgy^»  Is  given  by 

the  formulation  in  section  II-C  using  the  function  g(w).  The  real  and 
imaginary  parts  of  the  g(oj)  in  the  present  section  are  to  be  identified  as  the 
real  function  g(cj)  of  section  II-C.  The  equation  to  be  solved  is  equation  (60). 
The  choice  of  uncertainty  matrix  (i.e.  (V®))  for  the  data  values,  G(w^.,  is 
made  as  described  in  appendix  D.  The  £ matrix  and  the  V-vector  of  equation  (60) 
are  replaced  with  the  new  £ matrix  and  the  new  V-vector  of  equation  (79).  The 
derivative  terms  (3g^/3aj)  of  equation  (60)  are  discussed  in  another  section 
of  this  report. 
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Figure  20  shows  the  sequence  of  curves,  (r(w)J|^,  that  are  the  result  of 
the  reflection  coefficient  computations  made  using  each  (kj)  new  electron 
density  profile  obtained  from  the  iterative  scheme  as  applied  to  solutions 
of  equation  (60).  It  is  observed  in  the  figure  that,  although  the  sequence 
of  computed  point  r(cj^)  in  the  curve  [r(w)J^  is  completely  reversed  from  that 
of  the  data  points  R(wj),  the  iterative  procedure  causes  the  successive  curves 
(i.e.  Irioj)]^,  (i"(‘*j)]3  etc.)  to  undergo  a transformation  (i.e.  counterclockwise 
rotation)  that  leads  to  the  desired  match  between  computation  and  data. 

Note  that  this  choice  of  function  g(uj)  to  replace  r(co)  is  an  intuitive  one, 
and  that,  at  best,  probably  only  one  sense  of  rotation  (clockwise  vs.  counter- 
clockwise) will  avoid  the  singular  coefficient  matrix  of  equation  (60).  An 
example  has  been  found,  however,  with  the  use  of  simulated  data,  when  a correct 
choice  of  rotation  sense  led  to  convergence  in  the  sequence  of  solutions  to 
equation  (60). 

H.  The  Derivatives  (dr/dw)^  and  (dR/dw). 

There  is  a requirement  for  a method  to  define  the  derivatives  (dr/dw)^ 
and  (dR/dw)^  along  the  curves  r(w)  and  R(w)  when  values  are  given  only  at 
the  data  frequencies  These  given  values  are  r(^^)  and  R(t*J!.)  with 
i = l,...,m.  A way  to  satisfy  this  requirement  is  to  use  the  interpolation 
scheme  described  in  section  II-B-2.  The  derivative,  (dr/dw)^  ( (or  dR/dw)^] , 
is  defined  in  terms  of  a smooth  curve  drawn  through  the  points  r(w^)  (or  R(t*>^)I 
by  the  interpolation  scheme. 

The  smooth  curve  is  represented  in  terms  of  closely  spaced  points  p(<^|()> 
at  an  interval  in  frequency  of  Aw,  which  include  the  set  of  r(w^)*s 


More  precisely,  the  interpolation  procedure  assumes  that  values  of 
Wj^are  chosen  at  k = These  values  are  separated  by  the  incre- 

ment, Aw.  It  is  also  assumed  that  known  values  of  the  parameters  r(aj^) 
(or  R(u).]  exist  at  the  rvalues,  where  i=l,...n.  Note  that  K will  be 
greater  than  n.  In  the  set  of  there  will  exist  values  of  w such  that 
each  w^will  correspond  to  an 

In  section  II-B-2,  it  is  shown  that,  using  the  interpolation  scheme 
described,  the  variable  =p(W|^)may  be  written  as  a linear  combination  of 
the  values  of  the  r(co^.)'s  [or  the  R(co.)'s].  This  is  expressed  as  in 
equation  (32)  as: 


At  these  values  of  cj  where  an  is  equal  to  an  the  value  of  p{.o)y)  will 
be  equal  to  the  r(Gj^)  [or  R(oj^)]  value. 

The  derivative  at  one  interpolation  point,  is  defined  in  terms  of 
the  change  in  value  of  p(w)at  the  interpolation  points  preceding  and  follow- 
ing that  point.  That  is: 

a ' &CO 

where  Aw  is  the  interpolation  increment  in  frequency  and  the  p(W|^)'s  are  the 
interpolated  values. 
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It  Is  also  shown  In  section  II-B-2,  that  the  derivatives,  (dp(w)/dcj)|^ 
are  linear  combinations  of  the  values  of  r(a;^)[or  R(w^.)l-  That  Is: 
from  equation  (35): 


I 

1.  ->• 


In  particular,  at  the  data  frequencies,  , 


and  similarly: 


where 

for  values  of  "k''  where  corresponds  to  the  data  frequencies,  . 
That  Is,  since: 
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and 


Then,  only  those  rows  of  the  £ matrix,  which  correspond  to  those  rvalues 
where  u)|^=w^,are  substituted  into  equation  (96)  in  the  row  positions  to 
give  values  to  the  a matrix. 

For  a given  value,  equation  (92)  may  be  written: 


Note  also,  that  derivatives  of  (dr/d  w) ^ wi th  respect  to  the  electron 
density  profile  parameters,  aj,  are  of  the  form: 


L(dL) 


The  above  procedure  of  determining  values  for  (dr/dw)^  and  (dR/dw)^  allows 
for  the  phase  terms  of  equations  (84)  and  (86)  to  be  followed  along  the 
curves  r(w)  and  R(cj)  so  as  to  ensure  continuity  of  phase  in  selecting  the 
values  of  the  n^  parameters  in  the  functions  g(u^)and  G(co^). 
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I.  Following  of  the  Phase  Angle,  0.  Along  the  Curves.  rM  [or  RMl  , 
and  Derivation  of  the  Derivatives,  (dg(u)j )/daj) 

The  curve,  r(t>;)  [or  R(w)]  of  section  III-F  Is  defined  In  terms  of  the 
Interpolation  points  of  section  III-H  at  frequencies,  spaced  Au  apart. 
From  equation  (84): 


and 

where  4>|^is  the  principal  part  and  the  nj^'s  are  chosen  so  that: 

<t>i,  - 

are  taken  as  small  along  the  g(w)  [or  G(w)]  curve  to  ensure  continuous  changes 
In  phase. 

At  the  data  frequencies,  w.,  equation  (84)  Is: 


= Ir,  ( j L = - • -1  (loi^ 


a 


nol 


Now  taking  the  derivative  of  equation  (101)  with  respect  to  the  electron 


density  profile  parameters,  a.,  gives: 


(io3) 


Note  that  the  number  of  cycles,  n-,  are  not  relevant  to  the  derivative. 

At  the  data  frequencies,  (cj .),  equation  (97)  of  section  III-H  may  be 
substituted  into  equation  (101)  above  to  get  for  each  g^: 


‘J.  - / 


m 


h 


\ n- . • r- 


then,  using  the  following  identity: 


ci  U 


d In  i ‘■0  

c<«  — U da 


the  derivative  (dg./da.)  may  be  written: 

* w 


5 9.  _ 


1 


^ ^ t 


f- 


(KP-f) 


(io^ 


where  the  (a^^i)  values  are  given  in  equation  (94)  of  section  TII-H  and  the 
dr^/doj)  terms  are  given  in  equations  (70-73)  in  section  III-E. 
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IV.  SUMMARY  OF  THE  INVERSION  TECHNIQUE  WITH  EXAMPLE  USING  SIMULATED  DATA 
As  stated  in  the  introduction,  the  approach  taken  in  finding  a best-fit 
electron  density  profile  from  data  is  an  iterative  technique  which  optimizes 


the  trade-off  between  the  deviation  from  data,  represented  by  the  value 
"s",  and  curvature  in  the  wanted  electron  density  profile,  represented  by 
the  value  of  "c".  The  trade-off  is  given  by  the  condition: 

Q — |—  ^ S C>  /TTtmmum 


where  "X"  must  be  chosen  for  the  optimum  condition. 


The  actual  equation  to  be  solved  in  order  to  obtain  the  desired  electron 
density  profile  is,  from  equation  (60): 


The  terms  (Cj.j),  (3g^/3«j),  (V),  (g^  - Gj),  (wj)®  and  (Aaj)are  all 

described  in  section  III. 

Equation  (108)  may  be  written: 


^ 


R ^ ‘‘O)  ('l 
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where 


The  coefficient  matrix,  and  is  a function  of  the  trade-off 
parameter,  "X";  and  also  a function  of  the  computed  ionospheric 
reflection  coefficient  r(cj). 

A vector  which  is  a function  of  the  above  "X",  and  r.  It  is  also 
a function  of  the  ionospheric  reflection  coefficient  data, 

R(cj),  and  a function  of  the  previous  electron  density  profile, 

A vector  which  represents  the  modifications  to  be  added  to  the 
previous  electron  density  profile  for  that  step  of  the  iteration 
scheme. 

The  data  frequency,  i=l,,..m. 

Equation  (109)  is  solved  for  (Aaj)as: 


Recall  that  (a.°)  refers  to  the  previous  electron  density  profile  in  the 

J 

iteration  scheme  (e.g.  the  initial  profile  in  the  case  of  the  first  iteration 
step)  and  that  (Aa  .)  is  the  solution  to  each  iteration  of  equation  (108) 

(or  (110)1.  The  new  profile  generated  by  each  iteration  is: 

' -f-  (3^5 

where  each  term  is  a vector  of  length  j=l,...n,  with  n being  the  number  of 
segments  in  the  electron  density  profile. 

A full  wave  computation,  using  each  set  of  (aj)'s  as  the  input  profile, 
will  give  the  set  of  reflection  coefficients,  r(w^)  needed  in  the  inversion 
scheme. 


K = 


T = 


(i«j)  ' 


CO  . = 

1 
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An  example  of  the  inversion  procedure  is  illustrated  in  figures  21, 

(a  and  b)  and  22,  (a  and  b).  In  particular,  figure  21, a shows  a data  profile 
for  which  simulated  reflection  coefficient  data  was  obtained  by  a full -wave 
computation.  Also  shown  in  the  figure,  is  the  starting  profile  used  in  the 
inversion  technique.  The  positions  of  the  profile  "stops"  are  indicated  in 
the  figure. 

Figure  21, b gives  a comparison  between  the  simulated  R(cj^)data  values, 
as  obtained  from  the  data  profile,  and  the  computed  r(w^)  values  from  com- 
putations using  the  starting  profile  of  figure  21, a.  The  data  frequen- 
cies used  were  from  8 through  17  kHz  at  a frequency  increment  of  1 kHz. 

The  8 kHz  computed  value  is  identified  by  the  box  that  encloses  the  number  1 
while  the  17  kHz  value  is  identified  by  the  box  with  the  number  10.  In  this 
example  it  was  assumed  that  the  R(w^.)data  was  error  free.  The  figure 
illustrates  that  the  R(c<;.)  values  and  the  r(cjp  values  are  very  different  at 
this  initial  step  of  the  iteration. 

Figure  22, a shows  the  comparison  between  the  original  data  profile  and  the 
final  profile  obtained  from  the  inversion  scheme.  Note  in  the  figure  that 
the  top  of  the  final  profile  has  moved  away  from  the  original  "':top"  value, 
thus  indicating  that  information  exists  in  the  data  values  that  leads  to  this 
characteristic.  The  lowest  value  of  the  final  profile  remains  at  the  stop 
value  indicating  that  no  information  is  available  in  the  data  to  modify  this 
relationship. 

Figure  22, b illustrates  the  comparison  between  the  data  values  (R(w^) 
and  the  computed  values  (r(cj^)as  obtained  from  the  final  profile  of  figure 
22, a.  In  this  case  the  r{w)  curve  has  completely  reversed  in  sequence  from 
that  shown  in  figure  21, b.  This  was  accomplished  by  a counterclockwise 
rotation  of  the  computed  curve  r(cu)  as  the  iterative  sequence  progressed. 


Also,  the  spacing  between  the  various  r(co.)  values  has  been  modified  from 
that  of  figure  21, b so  that  the  fit  to  the  error-free  simulated  data,  R(w^)  , 
appears  to  be  very  good. 

It  must  be  realized  that  the  quality  of  fit,  shown  above,  between  the 
simulated  data,  ),  and  the  computed  values,  r(w^),  as  obtained  from  the 

inversion  procedure,  are  done  for  the  ideal  case  of  no  error  in  the  R(tj^) 
values.  Also,  it  was  assumed  that  the  electron-neutral  particle  collision 
frequency  was  known  exactly.  In  cases  where  the  input  parameters  are  not  as 
described  in  the  above  example,  the  results  cannot  be  expected  to  be  as  out- 
standing. Examples  of  other  cases,  in  which  gaussian  random  error  is  added 
to  the  simulated  data,  are  given  in  reference  8 which  also  describes  the 
computer  program  which  applies  the  inversion  scheme  described  in  this  report. 
The  quality  of  results  obtained  using  actual  experimental  measurements  is  still 
to  be  determined. 
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Figure  8.  Comparison  of  optimally  smoothed  and  undersmoothed  results  to 
data  with  added  gaussian  error. 


Comparison  of  optimally  smoothed  and  greatly  undersmoothed  results  to  data  with  added  gaussian  error. 


Figure  12.  Increments  associated  with  the  weighting  funclion  (da/dajl. 


ERROR  FREE  DATA 


Figure  13.  R-plane  represenlation  of  R (real) 
and  R (imaginary)  data  as  a function  of 
propagation  frequency,  uij. 


Figure  14.  Constraints  on  the  electron  density  profile. 


Figure  15.  r-plane  representation  of  r(uj)  (real) 
and  r(u))  (imaginary)  as  a function  of 
propagation  frequency. 


R,  r IrMli 


Figure  1 7.  Matclied  comparison  of  R(bJi)  reflection 
coefTicknl  data  and  r(cj)  compuUtions. 


Figure  16.  Comparison  of  R(u)i)  reflection 
coefficient  data  and  rfw)  compulations. 


R,  r (real) 


Figure  18.  Possible  sequence  of  r(u>)  iterations  to  obtain 
match  between  and  R(Wj)  values  in  which  a “cusp" 
is  allowed  to  form. 


Figure  l‘).  The  r(w)  curve  illustrating  the  transl'orniation  function,  glcjj)  = In 
(i.e..  spacing  and  phase). 


Figure  20.  Sequence  of  r(co)  iterations  to  obtain  match 
between  r(wi)  and  R(u)i)  values  using  the  g-function. 


t R,  r UdaO 

1 

i 

Figure  22.  Inversion  procedure,  final  iteration,  for  error-free 
simulated  data. 
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VII.  APPENDICES 

Appendix  A.  The  "C"  Matrix 

In  the  procedure  of  smoothing  out  error  in  data  of  the  form  y(x)  vs. 
X,  a term  identified  as  curvature  is  utilized. 

The  curvature  term  is  defined  as: 


Equation  (A-1)  refers  to  a continuous  curve  y(x).  However  in  practice,  y(x) 
can  only  be  specified  at  a finite  number  of  points,  Xj  (j=l,...n).  For 
simplicity  it  will  be  assumed  that  the  x.'s  are  equally  spaced  at  intervals, 

J 

Ax. 

At  the  position  x.  the  curvature  term  is  approximated  by: 


d- 

dx 


\ Av  A.V  //  ' 


This  term  is  Illustrated  in  figure  A-1. 

The  curvature  function,  "c",  is  then  given  by  the  expression: 


h=2 


where  n is  the  number  of  points,  Xj. 
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Figure  A-1.  Plot  of  y(x)  for  curvature. 


Note  that  the  expression: 

Cvy.l  -2  1)^ 


can  be  written  as: 


For  illustration,  let  n ■ 8 end  expend  the  expression  (A-4).  This  gives, 
from  equation  (A-3);  (5iX^*-C  =• 

( ‘K  vj KyJ - Uu)  - ¥ ( yj.) + a ( y \') 

+( H Yj+  (yj-y  Oj3  - y (y,y,)ta  ( yj,) 

+( Y/^  y ( y,)'+  (yj-  y Cy,  y,)  - ( \ y.) + 2 C y,  y.) 
■^(  Yr)'+  y(  yj+  (yj-y  6,YrJ>~Y(rj.')-t-a(y,yj 

+C  v(  yj+  (yJ  - Y Oriy  cyj,) + ^ ( Y.x) 

-YCyJ+YC  1I,T+  (Yj-YCYA)-Y(Y7y,)+i(Yj,) 

Rearranging  terms  gives 


XI 

xx 

+ 'lA-X/sV^XX-^XX'^  1A 

^ Y,yr‘Y7J,+  ^Y,Y,-YY,Y.4  XX 

xxr^txx^ (>xXr-nx^-  ix 

" Y.rYiY/^UrYU+  XY,  : 

+ Y,y^-y^);+y'i;r,-jy,y,  i 

XX- ^XX*  XX,  ! 

9 


In  matrix  notation  equation  (A-6}  may  be  written  as; 

(ax)'-c  = 

1 -i  I o o O o o\ 
-Z  S'  -y  \ o O O o ' 

I -y  i -V  I o c £> 

0 / -1  6 -f  I o 

o 0 \ i-  -^  \ c 
c o o I -y  i / 
c o o o I -V  f -t 
c o o c I -t  I 


\ 


or 


i4)0-c= 


where  J ' 
or  j' 
Also 


1 9 • • • n and 

1 ....  8 or 


Ca-^ 


CA-s) 


J “ lf»»*n 
j = 1,...8 


^ i X S Vft-tvr 


/s  a/7  ^ X S wcftrDi 


/S  a 2 X ^ CChn^n  vectc>r 
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Therefore  the  derivative  terms  of  "c"  may  be  written  as; 


where  J'  = 
and  j = 1,...8 

It  Is  of  Interest  to  note  that  the  second  derivative  of  "c"  with  respect 
to  y.  Is: 


Appendix  C.  The  Derivative  Terms  (dy/dx)j 

Given  a set  of  data  values  Y(X^),  (1  = at  each  point  X^.  It 

is  required  to  find  a set  of  interpolated  values  y(Xj).  (j  = l,...n),  where 
the  Xj  are  evenly  spaced  increments.  Ax,  of  x.  It  is  required  that  the  set 
of  y(Xj)  values  include  the  set  of  Y(X^)'s.  That  is  each  Y(X^)  must  also 
be  a y(Xj). 

From  equation  (29)  of  section  II-B-2,  the  relationship  between  the 
Y(X^)'s  and  y(Xj) 's  is: 


where 

y,  = '/(X), 

where  the  j'-th  row  of  is  the  j'-th  row  of  the  matrix  £ of 

equation  (22),  If  y(Xj)  does  not  correspond  to  an  Y(X^).  If  y(x^)  does 
correspond  to  an  Y(X^),  then  the  j'-th  row  of  Is  all  zeros  except 

for  a 1 for  the  main  diagonal  element,  Cj.j,.  Also  all  the  elements  of  the 
vector  Vj.  are  zero  except  that  if  y(Xj)  corresponds  to  an  Y(X^),  then 


Vji  > Y(X^).  Equation  (30)  of  section  II-B-2  is  an  example  of  these 
relationships. 


Now  rewrite  equation  (C-1)  as 


where 


and  where 


and  also 


so  (Xj'j)  chosen  to  be  such 


that  it  contains  no  Y/s. 
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Now  consider  equation  (C-2)  with  the  vector,  CY^)  omitted  from  both  sides 
of  the  equation.  This  is: 


L—. 

tlrA 


VX^CcV/  _ V-  - 


which  contains  no  Y^'s,  but  which  may  be  solved  for  bj^,  that  is: 


i-  K,,  - * c' A ■ * 

' J i Wij  J ^ 


Once  the  (bj^'s)  are  known,  the  y^'s  may  be  found  from  equation  (C-3). 


Cc'-g) 


This  equation  illustrates  that  the  interpolated  values,  y^,  may  be 


found  from  a linear  combination  of  the  data  values,  Y^.  That  is: 


- f' 


Now  consider  a way  of  defining  (.dy/d;(}j  in  terras  of  y^  values.  Let 


•5.  AX 


Since  from  equation  (C-9) 


or 


J.  AK 


This  equation  may  be  written  as: 


where 


Xu  X 


76 


Appendix  D.  The  Uncertainty  Matrix,  [l/o] 


Given  equation  (54)  of  the  main  text  as: 


/\i  o\ 

where  I Is  a diagonal  matrix  with  real  constant  elements.  Also 

is  a data  value,  r^  is  a computed  value  of  the  curve  r(c<^)  and 
AOj  = “(Zj)  -a°(zj). 

Consider  the  left  side  of  equation  (D-1).  This  may  be  written  as: 
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Also,  consider  the  following  products  of  equation  (D-2): 


I icin’  i 


At  this  point  equation  (D-2)  has  not  been  modified.  The  equation  will  be  modi 
fled,  however,  by  a change  applied  to  the  product  factors  of  (D-5)  where: 


This  means  that  in  the  function  (^g^V^r.)  of  (D-5),  everywhere  an  r^ 
occurs  replace  it  with  (i.e.,  a data  value). 

Now  using  equations  (D-5),  {D-6)  and  (D-7)  define: 


where  equation  (D-8)  and  (D-9)  contain  only  constants. 

Substitution  of  equations  (D-3)  through  (D-9)  into  equation  (D-2)  gives: 


Note  that  as  the  parameter  "X"  is  increased  in  steps,  the  values  of  the  r^'s  of 
equation  (D-2)  approach  the  values  of  the  data  values,  R. . That  is,  in  equation: 
(D-3)  and  (D-4): 


i 


79 


Approaches 


as  "X"  becomes  larger. 

Thus  the  following  product,  implicit  in  equation  (D-lQl, 


approach  the  unit  matrix  as  the  r. 's  approach  the  R^. 's.  Then,  equation  (D-10) 
will  approach  being  the  original  left  hand  side  of  equation  (D-1). 

Now  consider  the  right  hand  side  of  equation  (D-1).  This  may  be  written 
as: 
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Note  in  equation  (D-14)  that: 


— 1!!». 


And 


Also  consider  the  following  products  of  equation  (D-14): 


- // 


At  this  point,  equation  (D-14)  has  not  been  modified,  except  by  way  of  approxi- 
mation (D-16).  Equation  (D-14)  is  modified  by  the  change  applied  to  the  product 
factors  of  (D-17)  where: 


Again,  this  means  that  in  the  function  (3g^./3r^)  of  (D-17),  everywhere  an 
r^  occurs  replace  it  with  (i.e.,  a data  value). 

Now  using  equations  {D-17),  (D-18)  and  (0-19)  define: 


where  (D-20)  and  (D-21)  contain  only  constants. 
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Substitution  of  equations  (D-15)  through  (D-21)  into  equation  (D-14)  gives: 


As  discussed  previously  for  the  left  hand  side  of  equation  (D-1),  as  the  param- 
eter "X"  is  increased  in  steps,  the  values  of  the  r^. 's  of  the  right  hand  side  of 
equation  (D-1)  approach  the  values  of  the  data  values,  R.. 

That  is,  from  equation  (D-15)  and  (D-16) 


approaches: 


as  "X"  becomes  larger. 

Thus  the  following  products,  impli 


c ^ 

{-'-K 

> > 
'•  -> 

u 

i 

it  in  equation  (D-14): 


approach  the  unit  matrix  as  the  r^'s  approach  the  R^'s.  Then  equation  (D-22) 
will  approach  being  the  original  right  hand  side  of  equation  (D-1). 
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Appendix  E.  The  Derivatives,  (ag^'/dr.). 

Appendix  0 makes  no  assumptions  about  the  functional  relationships  between 
the  vectors  T and  "g  (and  hence  R and  G)  except  that  the  derivatives  (3g^'/3r^) 
must  be  definable. 


From  equation  (97)  of  the  main  text: 


Substituting  equation  (E-1)  into  equation  (83)  gives: 

’ c'  - m 

Taking  derivatives  with  respect  to  r^.  and  using  the  identity: 

<A  jin  i t-O)-  I 

cJx  ~ M tU 


(f  ^ 


Gives  from  equation  (D-2): 


^ 1 


3 <J  Ci'i. 


ill  - 


^ -'1  - o 

L = ' ' i-n 


0 


(e-.) 


t"rl. 


and  similarly  for  R and  G. 
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Wright-Patterson  AFB,  OH  45433 
OICY  ATTN  ETD  B.  L.  Ballard 

HQ  USAF/RD 

Washington,  DC  20330 
OICY  ATTN  RDQ 
Headquarters 

North  American  Air  Defense  Command 
1500  East  Boulder 
Colorado  Springs,  CO  80912 
OICY  AHN  Chief  Scientist 

Conmander  , 

Rome  Air  Development  Center,  AFSC 
Griffiss  AFB,  NY  13440 

OICY  ATTN  EMTLD  Doc.  Library 

Commander 

Rome  Air  Development  Center,  AFSC 
Hanscom  AFB,  MA  01731 

OICY  ATTN  EEP  John  Rasmussen 

SAMSO/W 

Norton  AFB,  CA  92409 
(Minuteman) 

OICY  AHN  NTfiL  LTC  Kennedy 

Commander  in  Chief 
Strategic  Air  Command 
Offutt  AFB,  NB  68113 
OICY  ATTN  NRT 

OICY  ATTN  XPFS  MAJ  Brian  G.  Stephan 
OICY  ATTN  DOK  Chief  Scientist 


Other  Government 


US  Energy  Research  and  Dev.  Admin. 

Department  of  Energy 
Albuquerque  Operations  Office 
P.O.  Box  5400 
Albuquerque,  NM  87115 

01 CY  ATTN  Doc  Con  for  D,  W.  Sherwood 


Department  of  Commerce 
National  Bureau  of  Standards 
Washington,  DC  20234 
01  CY  AHN  Arthur  Ernst 
01 CY  ATTN  Raymond  T.  Moore 


Department  of  Energy 
• Division  of  Headquarters  Services 
Library  Branch  G-043 
Washington,  DC  20545 

01 CY  ATTN  Doc  Con  for  Allen  Labowitz 

Division  of  Military  Application 
Department  of  Energy 
Washington,  DC  20545 

01 CY  ATTN  Doc  Con  for  Donald  I.  Gale 

University  of  California 
Lawrence  Livermore  Laboratory 
P.O.  Box  808 
Livermore,  CA  94550 

01 CY  ATTN  Glenn  C.  Werth  L-216 
01  CY  ATTN  Tech.  Info  Dept  L-3 
« 01  CY  AHN  Frederick  D.  Seward  L-46 


Department  of  Commerce 
Office  of  Tel econinuni cations 
Institute  for  TELCOM  Science 
Boulder,  CO  80302 

01  CY  AHN  William  F.  Utlaut 
01 CY  ATTN  L.  A.  Berry 
01 CY  ATTN  A.  Glenn  Jean 
01 CY  ATTN  D.  D.  Crombie 
01 CY  ATTN  J.  R.  Wait 

Department  of  Transportation 
Office  of  the  Secretary 
TAD-44.1,  Room  10402- B 
400  7th  Street,  SW 
Washington,  DC  20590 
01  CY  AHN  R.  L.  Lewis 
01  CY  AHN  R.  H.  Doherty 

Department  of  Defense  Contractors 


Los  Alamos  Scientific  Laboratory 
P.O.  Box  1663 


Los  Alamos,  NM  87545 

01  CY  ATTN  Doc  Con  for  T.  F. 

01  CY  AHN  Doc  Con  for  D.  R. 

01 CY  ATTN  Doc  Con  for  P.  W. 

OICY  ATTN  Doc  Con  for  J.  H. 


Taschek 

Westervelt 

Keaton 

Coon 


Sandia  Laboratories 
Livermore  Laboratory 
P.O.  Box  969 
Livermore,  CA  94550 

OICY  ATTN  Doc  Con  for  B.  E.  Murphey 
OICY  ATTN  Doc  Con  for  T.  B.  Cook  ORG  8000 


Aeronony  Corporation 
217  S.  Neil  Street 
Champaign,  IL  61820 
OICY  ATTN  S.  A.  Bowhill 

Aerospace  Corporation 

P.O.  Box  92957 

Los  Angeles,  CA  90009 

OICY  ATTN  Irving  M.  Garfunkel 

Analytical  Systems  Engineering  Corp. 
5 Old  Concord  Rd. 

Burlington,  MA  01803 
OICY  ATTN  Radio  Sciences 


Sandia  Laboratories 
P.O.  Box  5800 
Albuquerque,  NM  87115 

OICY  ATTN  Doc  Con  for  Space  Proj  Div 
OICY  ATTN  Doc  Con  for  A.  D.  Thornbrough 
ORG  1245 

OICY  ATTN  Ox  C n for  W.  C.  h^yra 

OICY  ATTN  Dx  Con  for  3141  Sandia  Rpt  Coll 


Boeing  Company,  The 
P.O.  Box  3707 
Seattle,  WA  98124 

OICY  ATTN  Glenn  A.  Hall 
OICY  ATTN  H.  R.  Millard 
OICY  ATTN  J.  F.  Kenney 
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University  of  California 
at  San  Diego 

Marine  Physical  Lab  of  the 
Scripps  Institute  of  Oceanography 
San  Diego.  CA  92132 

01 CY  ATTN  Henry  G-  Booker 


Geophysical  Institute 
University  of  Alaska 
Fairbanks,  AK  99701 
01 CY  ATTN  T.  N.  Davis 
01 CY  ATTN  Neal  Brown 
01 CY  Technical  Library 


ConY)uter  Sciences  Corporation 
P.O.  Box  530 
6565  Arlington  61 vd 
Falls  Church,  VA  22046 
01 CY  ATTN  D.  Blumberg 


GTE  Sylvania,  Inc. 
Electronics  Systems  GRP 
Eastern  Division 
77  A Street 
Needham,  NA  02194 

01  CY  AHN  Marshal  Cross 


University  of  Denver 
Colorado  Seminary 
Denver  Research  Institute 
P.O.  Box  10127 
Denver,  CO  80210 

OICY  ATTN  Donald  Dubbert 
01  CY  ATTN  Herbert  Rend 

Develco 

530  Logue  Avenue 
Mountain  View,  CA  94040 
OICY  ATTN  L.  H.  Rorden 

ESL,  Inc. 

495  Java  Drive 
SuMVvale.  CA  94066 

OICY  ATTN  James  Marshall 

General  Electric  Company 
Space  Division 
Valley  Forge  Space  Center 
Goddard  Blvd  King  of  Prussia 
P.O.  Box  8555 
Philadelphia.  PA  19101 
OICY  ATTN  M.  H.  Bortner 
Space  Science  Lab. 

General  Electric  Conqjany 
TEl^'Center  for  Advanced  Students 
816  State  Street 
P.O.  Drawer  QQ 
Santa  Barbara,  CA  93102 
OICY  ATTN  B.  Gwibill 
02CY  ATTN  DASIAC 
OICY  ATTN  Don  Chandler 
OICY  ATTN  Warren  S.  Knapp 


IIT  Research  Institute 
10  West  35th  Street 
Chicago,  IL  60616 

OICY  ATTN  Technical  Library 

University  of  Illinois 

Department  of  Electrical  Engineering 

Urbana,  IL  61803 

02CY  ATTN  Aeronony  Laboratory 

Johns  Hopkins  University 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel,  MD  20810 

OICY  ATTN  Document  Librarian 
OICY  ATTN  J.  Newland 
OICY  P.  T.  Komiske 

Lockheed  Missiles  & Space  Co,  Inc. 
3251  Hanover  Street 
Palo  Alto,  CA  94304 
OICY  ATTN  E.  E.  Gaines 
OICY  ATTN  W.  L.  Imhof  D/52-12 
OICY  ATTN  J.  B.  Reagan  D/52-12 
OICY  ATTN  R.  G.  Johnson  D/52-12 

Lowell  Research  Foundation,  Univ.  of 
450  Aiken  Street 
Lowell,  MA  01854 
OICY  ATTN  Dr.  Bibl 

M.I.T.  Lincoln  Laboratory 
P.O.  Box  73 
Lexington,  MA  02173 
OICY  ATTN  Dave  White 
OICY  ATTN  J.  H.  Pannell  L-246 
OICY  ATTN  R.  Enticknapp 
01  a ATTN  0.  M.  Towle 
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Mission  Research  Coporation 
735  State  Street 
Santa  Barbara,  CA  93101 
01 CY  ATTN  R.  Hendrick 
01 CY  AHN  F.  Fajen 
01  CY  AHN  M.  Scheibe 
01 CY  ATTN  J.  Gilbert 
01  CY  ATTN  C.  L.  Longtnire 

Mitre  Corporation 
P.O.  Box  208 
Bedford,  MA  01730 
01  CY  ATTN  G.  Harding 

Pacific-Sierra  Research  Corp. 

1456  Cloverfield  Blvd. 

Santa  Monica,  CA  90404 

01  CY  AHN  E.  C.  Field,  Jr. 

Pennsylvania  State  University 
Ionosphere  Research  Laboratory 
318  Electrical  Engineering  East 
University  Park,  PA  16802 

02CY  ATTN  Ionospheric  Rsch  Lab. 

R&D  Associates 
P.O.  Box  9695 
Marina  Del  Rey,  CA  90291 
01 CY  ATTN  Forrest  Gilmore 
01 CY  ATTN  William  J.  Karzas 
01 CY  ATTN  Phyllis  Grei finger 
01 CY  ATTN  Carl  Grei finger 
01  CY  AHN  H.  A.  Ory 
01 CY  ATTN  Bryan  Gabbard 
01  CY  AHN  R.  P.  Turco 


Rand  Corporation 
1700  Main  Street 
Santa  Monica,  CA  90406 

02CY  ATTN  Technical  Library 
01 CY  ATTN  Cullen  Crain 

SRI  International 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 
01  CY  ATTN  E.  T.  Pierce 
01 CY  ATTN  Donald  Neil  son 
01 CY  ATTN  George  Carpenter 
01  CY  AHN  W.  G.  Chestnut 
01  CY  AHN  J.  R.  Peterson 
01 CY  ATTN  Gary  Price 

Stanford  University 
Radio  Science  Laboratory 
Stanford,  CA  94305 

01 CY  ATTN  R.  A.  Helliwell 
01 CY  ATTN  R.  Fraser- Smith 
01 CY  ATTN  J.  Katsufrakis 

TRW  Defense  & Space  Sys.  Group 
One  Space  Park 
Redondo  Beach,  CA  90278 
01 CY  ATTN  Saul  Altschuler 
01 CY  ATTN  Dianna  Dee 


